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SECTION  1 
INTRODUCTION 


Nuclear-induced  lightning  (NIL)  has  been  observed  on  several  occa¬ 
sions  in  atmospheric  tests  in  the  1950s  and  1960s.  Perhaps  the  best  docu¬ 
mented  shot  was  Shot  Mike  of  Operation  Ivy  on  Eniwetok  Atol  on  October  31, 
1952  (Ref.  1).  Shot  Mike  was  a  10.4  MT  detonation  which  triggered  at  least 
5  lightning  flashes  which  grew  from  the  region  of  a  construction  camp  on  one 
side  and  from  water  on  the  other  »1  km  from  ground  zero,  to  heights  in  ex¬ 
cess  of  600  m  where  their  further  development  was  obscured  by  cloud  cover. 
One  frame  from  a  movie  is  shown  in  figure  1.  These  were  first  noted  at 
«1  ms  after  zero  time  and  persisted  for  «75  ms.  Calculations  would  suggest 
that  the  electric  fields  in  the  vicinity  of  these  lightning  strikes  were 
about  3  x  104  V/m  and  that^  the  dose  in  this  vicinity  was  =>3  x  107  rads/s 
(Ref.  2).  Propagating  strikes  were  observed  with  high  speed  photography  and 
from  these  records  one  may  infer  a  propagation  velocity  of  «1  to  2  x 
107  cm/s. 

Nuclear-induced  lightning  appears  to  differ  from  natural  lightning 
in  some  respects.  First  NIL  appears  to  grow  at  a  nearly  constant  rate  per¬ 
sisting  for  10's  to  100's  of  milliseconds,  whereas  natural  lightning's 
growth  is  discontinuous,  exhibiting  peak  currents  of  »104  A  developing  in 
times  of  »10-7  s  (I  “lO11  A/s)  with  a  fall  time  for  the  first  stroke  of 
®50  us.  Second  NIL  appears  to  propagate  at  a  slower  velocity  and  in  lower 
electric  fields.  Third  there  is  no  evidence  of  a  stepped  leader  which  is 
normally  associated  with  natural  lightning.  Although  nuclear  lightning 
transports  lower  peak  currents  than  natural  lightning  it  persists  so  long 
that  it  delivers  more  total  energy  to  the  lightning  rod.  (The  continuing 
current  in  natural  lightning  may  persist  for  seconds  at  very  low  levels.) 
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The  significance  of  NIL  is  that  it  may  represent  an  electromag¬ 
netic  threat  to  systems  that  would  otherwise  survive.  For  example,  NIL 
might  couple  to  antennas  or  the  power  distribution  systems  of  blast-hardened 
shelters. 


While  NIL  has  long  been  known,  it  was  only  recently  that  a  model 
was  created  to  explain  its  characteristics  (Ref.  2,3).  In  this  model,  elec¬ 
trostatic  fields  are  produced  by  photoCompton  currents  which  are  driven  from 
the  detonation  point  by  ionizing  radiation  and  by  the  return  current  through 
the  comparatively  conducting  earth  or  sea  surface.  While  these  electric 
fields  are  not  sufficient  to  produce  ordinary  electrostatic  air  breakdown 
(»3  x  106  V/m),  they  are  sufficient  to  produce  large  current  densities  in 
air  which  has  been  rendered  conductive  by  ionizing  radiation.  If  one  imag¬ 
ines  an  initiation  point  (such  as  the  antenna  or  structural  towers  observed 
in  Shot  Mike),  one  may  calculate  a  radius  within  which  the  electric  fields 
achieved  an  avalanche  value  ("3  x  106  V/m).  Whereas  the  electric  fields  are 
fixed  within  this  radii,  the  current  density  increases  as  r-2 ,  so  that 
within  a  smaller  region  of  radius  "a"  the  inwardly  flowing  power  density  is 
sufficient  to  (adiabatical ly)  heat  the  enclosed  air  to  a  temperature  which 
will  support  thermal  ionization  (»1  eV  or  11,200°C).  Thermally  ionized  air 
produces  a  small  conductive  bubble  which  subsequently  runs  into  a  filament 
of  nearly  constant  radius  which  burns  its  way  along  at  a  rate  determined  by 
the  rate  air  is  thermally  ionized  at  the  tip  of  the  discharge  (Ref.  2,3). 
This  model  is  schematically  represented  by  Figure  2. 

When  the  present  investigation  was  undertaken,  there  were  still 
significant  uncertainties  concerning  the  validity  of  this  model.  We  judged 
that  if  the  model  were  valid  and  could  be  extended  to  laboratory  simulator 
conditions  (i.e.,  AURORA)  by  using  the  scaling  relations  given  by  the  model, 
that  the  credibility  of  the  model  would  be  substantially  enhanced.  Conse¬ 
quently  a  laboratory  experiment  was  devised  (Ref.  4).  Inasmuch  as  no  labo¬ 
ratory  simulator  can  provide  the  dose  rate  associated  with  an  atmospheric 
test,  over  the  dimensions  and  duration  of  an  atmospheric  test,  considerable 
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PROPAGATION  VELOCITY 
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For  r  >  R  (Ohmic) 

J  =oE  =  KDE  E<E0 

For  r  >  R  (Avalanche) 

E  =  Ea  J=oEa(R/r)2 
P  =  JEQ 

Inside  r  ~  a  (Thermal  Ionization) 

v  .  Pa  .  KD(EQR)2 
11  ~Qa 

Experimentally  EQR  »Vappl 

D  measured 
a  estimated 

Qax  =«  /<7V2dt,  a  =  KD 
Qax  ss/ivdt 


Figure  2.  NIL  model . 


scaling  was  required.  Reference  4  suggested  that  applying  *  100,000  volts 
across  a  spark  gap  in  air  exposed  to  ^lO11  rads/s  would  result  in  an  ava¬ 
lanche  region  with  a  radius  of  approximately  1.6  cm.  The  predicted  propaga¬ 
tion  velocity  depends  on  the  instantaneous  power  density,  assuming  a  propa¬ 
gation  velocity  of  «107  cm/s  (from  Reference  1)  and  an  AURORA  pulse  width  of 
«120  ns,  the  NIL  is  expected  to  be  no  more  than  »1  cm  long. 

The  AURORA  experiment  was  designed  to  duplicate  the  conditions 
thought  to  exist  at  the  tip  of  a  propagating  nuclear  lightning.  The  AURORA 
simulator  provided  a  radiation  pulse  to  create  the  conductive  air.  The 
PULSAR  model  5Q  (100  kV)  pulser  was  syncronized  to  the  AURORA  pulse  to  pro¬ 
vide  the  electric  field.  In  a  three  day  test  period,  the  experiment  was  set 
up,  executed  and  torn  down.  At  least  five  events  were  observed  that  exhib¬ 
ited  the  charactersitic  of  a  NIL.  Specifically,  that  they  originated  at  a 
point  of  high  field  concentration,  propagated  at  a  velocity  of  «2  x  107  cm/s, 
and  persisted  only  while  the  radiation  pulse  provided  conductive  air. 

Nuclear  induced  lightning  can  only  occur  during  the  radiation 
pulse.  Any  streamer  that  draws  current  before  or  after  the  radiation  pulse 
is  not  NIL. 

A  description  of  the  experiment  is  presented  in  Section  2,  Experi¬ 
mental  Design.  The  results  are  presented  in  Section  3.  They  are  divided 
into  two  parts:  diagnostics  and  nuclear  lightning  results. 


SECTION  2 
EXPERIMENT  DESIGN 


The  objective  of  this  experiment  is  to  characterize  the  behavior 
of  a  lightning  stroke  in  dosed  (conductive)  air.  In  this  experiment  we 
created  conditions  indentical  to  those  postulated  to  exist  at  the  tip  of  a 
NIL  streamer  over  (limited)  dimensions  determined  by  the  capabilities  of  the 
radiation  source  (dose  rate  and  pulse  width)  and  electrical  pulser  (voltage 
and  stored  energy).  The  electrical  pulser  was  carefully  synchronized  to  the 
AURORA  photon  pulse.  The  electrical  conditions  were  measured  with  voltage 
probes  and  current  sensors.  Radiation  conditions  were  determined  from  TLD 
dosimetry  (dose)  and  transient  "pulse  shape"  by  dose  rate  sensors.  The  NIL 
event  was  sensed  electrically  and  optically.  An  optical  system  (dubbed  a 
macro-scope)  was  used  to  measure  the  discharge  propagation  characteristics. 
This  system  had  5  fast  detectors  each  of  which  viewed  a  2.5  mm  wide  strip  of 
the  object.  Open  shutter  photographs  were  obtained.  In  this  experiment  we 
created  a  radiation  induced  electrical  streamer.  The  key  characteristics  to 
be  investigated  include: 

a)  The  propagation  velocity  of  the  NIL; 

b)  The  current  carried  by  the  NIL; 

c)  The  applied  field 

d)  The  boundary  layer  potential  developed  at  the  cathode; 

e)  The  dependence  of  a),  b),  and  c)  on  the  dose  rate; 
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2.1 


TEST  CHAMBER  DESIGN 


The  test  chamber  geometry  was  designed  to  produce  fields  suffi¬ 
cient  to  produce  an  avalance  region  large  enough  for  a  streamer  of  measur¬ 
able  length  to  occur.  A  Pulsar  Pulsepak  50Q  (100  kV,  200  J)  pulser  was  used 
to  create  the  required  fields.  The  50Q  was  triggered  by  the  Pulsar  10A 
(10  kV)  trigger  generator.  The  geometry  is  limited  to  one  where  100  kV  is 
sufficient  to  create  the  avalanche  field  and  200  J  is  sufficient  to  provide 
current  during  the  120  ns  AURORA  pulse.  The  geometry  in  Figure  3  meets 
these  requirements.  The  hemispherical  cathode  tends  to  concentrate  the  cur¬ 
rent  flowing  to  the  pin  and  reduce  the  current  flowing  to  the  ground  plane. 
The  electric  field  at  the  tip  of  the  pin  could  be  slightly  adjusted  by 
changing  the  distance  between  the  hemisphere  and  the  ground  plane.  The 
radius  of  the  tip  of  the  pin  is  1.2  x  10-2  cm  (0.005").  The  test  chamber 
was  a  plexiglass  box.  The  chamber  was  mounted  on  a  wooden  frame  to  place  it 
at  the  height  of  the  AURORA  radiation  "hot  spot".  A  photograph  of  the  cham¬ 
ber  mounted  on  the  structure  at  AURORA  is  shown  in  Figure  4. 

2.2  ELECTRICAL  MEASUREMENTS 

In  order  to  electrically  characterize  this  laboratory  simulation 
of  NIL  we  equipped  the  experimental  test  chamber  with  a  number  of  voltage 
probes  and  current  sensors.  This  measuring  apparatus  is  schematically  indi¬ 
cated  in  Figure  5.  Voltage  from  a  100  kV  Pulsar  model  50Q  pulser  was 
applied  between  a  0.025  cm  (0.010")  diameter  wire  anode  and  a  hemispherical 
25  cm  diameter  stainless  steel  cathode.  The  voltage  applied  to  the  cathode 
(which  was  negative  with  respect  to  ground)  was  measured  with  two  voltage 
probes.  The  first  voltage  probe  was  comprised  of  a  28  kq  resistor  string, 
made  from  twenty  1.4  kq  resistors  and  a  1.8  q  measurement  resistor.  The 
response  time  of  this  probe  is  “28  ns  (Ref.  5).  A  second  voltage  probe  with 
greater  sensitivity  and  faster  response  time  (“10  ns)  was  comprised  of  a 
5.4  kq  resistor  string  composed  of  forty  135  q  resistors  and  a  1.8  q 


Figure  3.  Test  chamber  geometry. 


U  S£ 


current  sensors 


measurement  resistor.  Because  a  non-uniform  voltage  drop  can  occur  at  the 
cathode  (sometimes  called  cathode  fall)  a  grid  Langmuir  probe  was  installed 
»2  cm  away  from  the  cathode.  The  Langmuir  probe  voltage  was  monitored  with 
another  28  kft,  1.8  ft  voltage  divider.  The  Langmuir  probe  is  described  in 
Appendix  1. 

Radiation  induced  photo-currents  in  the  probe  and  resistor  string 
were  considered  to  determine  the  probable  signal  to  noise.  Photocurrents 
generated  were  expected  to  be  intermediate  between  the  Compton  current  from 
the  probe  area  («10-10  Coul /rad-cm2)  and  the  replacement  current  in  the  co¬ 
axial  cable  («10-13  A-s/rad-m).  The  overall  length  of  the  voltage  probe  and 
the  associated  cabling  was  «  1  m.  The  peak  dose  rate  at  the  probe  was 
“101 1  rad/s.  Thus,  the  radiation-induced  currents  are  intermediate  between 
10  mA  and  1  A.  The  anticipated  signal  was  generated  by  «  3  A  across  the 
sensing  resistor,  which  was  larger  than  the  expected  radiation-induced  cur¬ 
rents. 

To  measure  currents  to  the  pin  anode  or  annular  ground  plane  we 
measured  voltages  across  low  resistance  resistors,  as  shown  in  Figure  6.  To 
be  sure  this  measurement  is  correct  one  must  assure  oneself  that  the  shunt 
resistance  of  radiation  dosed  air  is  large  by  comparison.  The  resistance 
between  the  large  annular  plate  and  the  ground  plate  due  to  the  conducting 
air  can  be  calculated  from. 


where  d  is  the  plate  separation,  A  is  the  plate  area  and  a  is  the  conducti¬ 
vity  of  the  dosed  air.  o  is  calculated  from. 


a  =  K  y  e  u/a 
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where  K  is  2.2  x  109  cm-3rad_l,  u/a  is  2  x  10-5  cm2/V  for  damp  air,  e  is  the 
electron  charge  (1.6  x  10-19  Coni)  and  y  is  the  gamma  dose  rate  (rad/s). 

For  a  gamma  dose  rate  of  anything  less  than  1  x  1011  rads/s  (maximum  ob¬ 
served  dose)  and  two  1  m  diameter  plates  separated  by  2  cm  ,  the  resistance 
between  the  plates  is  about  0.4  ft.  The  measuring  resistor  between  the 
plates  was  set  at  0.018  i i  so  that  more  than  95%  of  the  current  flowed 
through  the  resistors.  The  responsivity  of  the  sensor  is  55.6  A/V.  Since 
currents  -30  A  were  expected,  500  mV  signals  should  be  observed. 

The  resistance  between  the  smaller  pin  plate  and  the  ground  plate 
due  to  the  conducting  air  was  “100  ft.  The  pin  was  expected  to  draw  -100  to 
1,000  A.  To  generate  signals  between  1  and  10  volts,  the  resistance  between 
the  pin  plate  and  the  ground  plate  was  0.018  ft.  The  resistance  between  the 
plates  in  both  cases  was  accomplished  by  connecting  100  1.8  ft  resistors  in 
parallel.  The  voltage  across  the  plates  was  fed  via  RG-58  cables  to  the 
data  room  and  recorded  on  oscilloscopes. 

2.3  OPTICAL  MEASUREMENTS 

The  propagation  velocity  of  the  NIL  was  measured  optically.  The 
optical  system  can  be  divided  into  two  categories:  a)  light  collection  and 
delivery  to  the  data  room;  b)  image  processing  and  recording. 

The  light  was  collected  with  two  parabolic  mirrors  placed  at  posi¬ 
tions  around  the  NIL  90°  apart  as  shown  in  Figure  7.  The  mirrors  were 
mounted  on  mirror  mounts  to  allow  for  angular  adjustments.  The  mirror 
mounts  were  mounted  on  translation  stages  to  allow  for  movement  along  the 
optical  axis  to  focus  the  mirrors.  The  tip  of  the  pin  in  test  chamber  was 
located  1.8  m  above  the  floor  (the  height  of  the  AURORA  "hot  spot"),  with 
the  test  chamber  placed  as  close  to  AURORA  as  possible.  The  two  beams  were 
sent  via  routing  mrirors  (mounted  on  mirror  mounts  and  tripods)  through  the 
optical  port  into  the  data  room.  The  beam  routing  system  is  represented  in 


Plexiglass 


Transfer 


Mirror 


Object 


AURORA 


f/6  Collector 


Wood  Frame 


Figure  7.  Test  chamber  and  parabolic  col  1 ection/coll imation 
mirrors.  Light  from  the  lightning  falls  on  a  pair 
of  30  cm  diameter  parabolic  mirrors  which  collimate 
the  light  and  send  it  back  through  the  test  chamber 
to  45°  flat  transfer  mirrors  which  deflect  it  away 
from  AURORA  toward  the  imaging  optics. 


Figure  8.  The  optical  port  was  located  «1.8  m  off  the  floor.  The  beams 
were  reflected  from  the  optical  port  to  the  image  processing  and  recording 
system  in  the  data  room.  Seven  mirrors  were  required  to  collect  and  trans¬ 
fer  the  image  of  the  NIL  to  the  data  acquisition  room.  The  seven  mirror 
surfaces  are  thought  to  have  a  reflectivity  of  «  80%  each,  resulting  with  an 
overall  transmission  of  21%. 

The  image  of  the  NIL  was  recorded  both  photographically  and  elec¬ 
tro-optical  ly.  The  optical  components  of  the  imaging  system  are  shown  in 
Figure  9.  Five  avalanche  photodiodes  (APDs)  were  placed  in  a  horizontal 
line  1  cm  apart.  A  1  cm  diameter  lens  with  a  focal  length  of  0.8  cm  was 
placed  in  front  of  each  photodiode.  The  output  of  each  APD  was  fed  into  an 
amplifier  with  a  gain  of  13  to  15  dB  and  an  input  impedance  of  50  a.  These 
were  the  receivers  of  the  DNA/HDL  high  speed  fiber  optic  links.  The  output 
of  the  amplifiers  was  fed  via  RG58  cables  to  the  oscilloscopes.  The  rise 
time  of  the  APD-amplfier  combination  was  »1  ns.  The  published  spectral  re¬ 
sponse  of  the  APDs  is  shown  in  Figure  10. 

In  Figure  11  the  optical  components  of  one  beam  are  shown.  To 
evaluate  the  magnification  and  the  power  losses  in  the  system,  the  system 
was  evaluated  in  the  direction  of  light  travel  (from  left  to  right).  Light 
from  the  (1  cm  long)  source  entered  reflector  A  (represented  in  Figure  11  by 
a  lens)  set  a  distance  f/\  (192  cm)  from  the  source.  The  energy  from  the 
source  radiates  isotropically.  The  optical  efficiency  of  this  reflector  is 
the  ratio  of  light  captured  by  a  mirror  of  radius,  r,  (Ax  =  irr2)  to  the  area 
of  a  sphere  whose  radius  is  the  mirrors  focal  length  (A2  =  ^(fA)2,  (r  « 

15  cm,  fA  *  192  cm).  This  is  0.15%.  Lens  C,  placed  30  m  to  the  right  of 
reflector  A,  will  create  an  image  100  cm  to  its  right.  The  beam  between 
optics  A  and  C  is  collimated.  The  image  (L)  is  located  at  the  focus  of  lens 
C  by  the  thin  (lens  formula), 

1  +  1  =1 
T  TJ  f 
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Figure  8.  Overhead  view  of  the  optics  system  (macroscope) 
designed  to  form  unity  gain  image  of  the  nuclear 
lightning  in  the  data  acquisition  room  about 
70  feet  away. 


Parabolic  Mirror  Image 

Mirror  A  c  Plane 


CO 


03  00 
+->  S- 
C  O 

o  +~> 
rsi  U 

•f-  QJ 

S- 

O  OJ 

n=  o 


r^\ 


*  * 

'0>, 

i  i 


Figure  9.  Optical  system  components. 
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where  the  object  distance  (0)  is  finite  (°°)  for  collimated  light.  The  mag¬ 
nification  at  this  image  is  given  by 

fc 

M  =  —  =  0.52 
fA 

where  fc  and  f^  are  the  focal  lengths  of  lens  c  and  mirror  A.  The  magnifi¬ 
cation  was  set  to  unity  by  adjusting  the  distance  from  the  object  to  collec 
tor  A  and  the  distance  from  the  lens  to  the  focal  plane.  The  size  of  image 
at  the  image  plane  is  then  »1  cm.  The  beam  is  22  cm  in  diameter  as  it 
passes  through  lens  C,  which  is  4.4  cm  in  diameter.  From  the  ratio  of  the 
areas  of  the  beam  and  the  lens,  only  10%  of  the  energy  caught  by  the  para¬ 
bolic  reflector  goes  through  lens  C. 

The  reader  should  notice  that  all  the  optics  used  in  this  investi 
gation  were  “off  the  shelf".  This  resulted  in  a  less  than  optimally  effi¬ 
cient  system.  Should  more  efficiency  be  required  in  further  investigations 
we  are  confident  that  at  least  a  decade  more  collection  efficiency  is  fea¬ 
sible. 

After  passing  through  image  plane  I,  the  beam  enters  cylindrical 
lenses  Bx  and  By.  The  magnification  performed  by  these  lenses  is  4.0  in 
the  "X"  dimension  and  0.25  in  the  "Y"  dimension.  The  purpose  of  the  cylin¬ 
drical  lenses  is  to  give  each  detector  a  band  shaped  field-of-view,  narrow 
in  the  arc  propagation  direction  and  wide  perpendicular  to  the  propagation 
direction.  Since  the  two  cylindrical  lenses  are  orthogonal,  the  image- 
object  relationship  can  be  treated  independently  for  each  lens.  Assuming  a 
source  which  is  1  cm  in  both  dimensions,  the  image  in  plane  II  will  be  4  cm 
in  the  X  dimension  and  0.25  cm  in  the  Y  dimension.  Of  course,  the  NIL  was 
essentially  a  line  -  so  this  represents  the  maximum  excursion  of  the  image 
in  each  direction.  All  of  the  light  passing  thorugh  lens  C  will  also  find 
its  way  through  Bx  and  By. 


The  diameter  of  the  detectors  in  image  plane  II  is  0.8  mm.  As 
mentioned  earlier,  8  mm  focal  length,  1  cm  diameter  lenses  were  placed  in 
front  of  each  detector.  In  the  "Y"  dimension,  the  beam  is  4  cm  wide.  With 
5  detectors,  the  exposed  "area"  is  4  cm  (4  detectors).  Therefore,  in  the 
"X"  dimension,  «  25%  of  the  energy  falls  on  a  detector.  About  4%  of  the 
power  is  lost  by  reflection  at  each  surface  of  lenses  C,  Bx>  By  and  the 
1  cm  lenses  at  the  photodiodes,  resulting  with  a  loss  of  28%.  The  spectral 
transmission  of  the  lenses  is  shown  in  Figure  12.  The  transmission  of  each 
lens  is  »  90%,  resulting  with  a  transmission  of  66%  for  all  4  lenses.  We 
estimate  the  total  system  efficiency  is: 

percent  of  NIL  power  sent  to  each  detector  = 


(0.15%)  x 

(21%)  x 

(72%) 

(66%) 
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capture 

losses 
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(4.4%) 
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=  1.65  x  10-4% 
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Figure  12.  Lens  Spectral  Transmission  Curve. 


The  optical  systems  sensitivity  is  estimated  here  to  be  1.65  pW/W. 
The  responsivity  of  the  APOs  shown  in  Figure  10  indicates  an  average  respon- 
sivity  in  the  visible  range  of  40  A/W.  Combined  with  the  amplifier  gain  the 
APD  response  is  9.3  kV/watt  (watt  on^  the  detector).  The  overall  optical 
system  response  (from  source  power  to  APD-ampl i f ier  output  voltage)  is  thus 
estimated  to  be  “15  mV/W. 

Polaroid  film  was  a  placed  at  image  plane  I  of  beam  II.  The  spec¬ 
tral  sensitivity  (cm2/erg)  required  of  the  film  can  be  calculated  from: 


Pi 


where  D  is  the  area  of  the  image  on  the  film  (cm2),  i  is  the  pulse  length 
(seconds)  and  P  is  the  radiated  power  at  the  film  (erg/s).  For  a  1  cm  x 
0.02  cm  image,  D  is  2  x  10~2  cm2.  The  AURORA  pulse  is  120  ns  FWHM,  so  t  is 
“120  ns.  For  10  Watts  (10y  erg/s)  from  the  object  in  this  spectral  region 
the  power  at  the  film  (P)  is  10B  erg/sec  x  6.4  x  10~b  (system  transmi ssion) 
“640  ergs/s.  The  required  spectral  sensitivity  is  then  260  cmz/erg.  The 
spectral  sensitivity  of  Polaroid  Type  612,  ASA  20,000,  is  “1,000  cm2 /erg. 
This  is  substantially  more  sensitive  than  Type  47,  ASA  3,000  film  which  has 
a  spectral  sensitivity  of  « 100  cm2/erg.  Type  612  film  was  used.  This  cal¬ 
culation  suggests  Type  612  was  more  than  sufficiently  sensitive,  whereas 
Type  47  would  have  been  marginally  adequate. 

Variable  apertures  were  placed  in  the  system  to  limit  the  inten¬ 
sity  of  the  beam,  if  required.  All  lenses,  apertures  and  detectors  were 
mounted  on  X-Y  translation  stages  to  facilitate  centering  and  adjusting  the 
beam.  All  of  the  optics  were  mounted  on  aluminum  rails.  The  rails  were 
mounted  on  1"  thick  plywood  which  was  placed  on  a  sturdy  table.  The  Pola¬ 
roid  film  was  held  in  an  oscilloscope  camera  back.  The  metal  slide  was  used 


as  a  manual  shutter.  Optical  alignment  was  performed  by  placing  a  600  watt 
tungsten  iodide  bulb  at  the  arc  location  and  adjusting  each  element,  start¬ 
ing  with  those  nearest  to  the  bulb. 

2.4  OPTICAL  SYSTEM  CALIBRATION 

The  optical  system  was  calibrated  by  using  the  600  watt  tungsten 
light  bulb  with  a  1  cm  long  filament  as  a  source.  The  light  from  the  light 
bulb  passed  through  the  optical  system  in  precisely  the  same  manner  as  the 
light  from  the  NIL  (with  the  exception  of  a  50%  transmission  neutral  density 
filter  in  the  optical  path).  Since  the  APD's  amplifier  is  AC  coupled,  the 
calibration  of  the  system  using  a  continuous  light  source  must  be  performed 
without  using  the  amplifiers.  To  do  this,  the  current  through  the  50  n 
input  impedance  of  the  amplifiers  was  measured.  The  1  cm  long  filament 
source  was  projected  onto  the  lenses  in  front  of  four  photodiodes.  Each 
photodiode  received  25%  of  the  light  from  the  filament.  The  tungsten  fila¬ 
ment  operated  at  about  3400  K.  It  has  been  estimated  that  10%  of  the  lamp 
power  is  radiated  in  the  visible  spectral  region  (Ref.  6).  From  this  esti¬ 
mate  the  600  W  bulb  used  for  this  calibration  will  radiate  60  W  of  visible 
light  (the  remainder  is  mostly  IR).  The  sensitivity  of  the  system  can  be 
calculated  from  the  current  amplifier  gain,  bulb  power  and  neutral  density 
filter  losses: 

(Photodiode  current)x(Ampl ifier  input  impedance)x(Ampl i fier  voltage  gain) 

(0.50  N.D.  Filter  transmission)  x  (60  watts)  x  (0.25) 

sensitivity  (V/watt) 

The  results  of  this  measurement  are  shown  in  Table  1.  The  average  sensiti¬ 
vity  of  the  photodiodes  is  14.8  ±  3  mV /watt.  The  uncertainty  principally 
reflects  uncertainties  in  the  APD  and  amplifier  gains.  In  this  test  we  did 
not  independently  check  the  output  of  the  light  bulb  or  the  accuracy  of  the 


neutral  density  filter.  The  National  Bureau  of  Standards  provides  suitable 
services  for  a  price.  The  spectrum  of  the  arc  was  predicted  (Ref.  7)  to  be 
different  than  that  of  a  gray  body  radiation  (i.e.,  the  light  bulb).  The 
spectra  of  the  filament  reference  source,  and  the  NIL,  and  the  spectral  re¬ 
sponse  of  the  optical  system  must  be  considered  when  calculating  the  total 
radiated  power  from  the  NIL.  Care  must  also  be  taken  to  assure  that  the 
lenses  in  front  of  the  photodiodes  are  completely  filled.  This  is  assured 
if  signals  are  observed  on  the  two  adjacent  photodiodes. 


Table  1.  Optical  system  calibration  data. 

(600  W,  3400°K  tungsten  filament  source) 


APD 

Current 

Sensitivity 

1 

50  uA 

10.4  ±  2  mV/w 

2 

16.6  ±  3  mV/w 

3 

95  uA 

19.7  ±  3  mV/w 

4 

60  uA 

12.4  ±  2  mV/w 

2.5  FACILITIES,  TIMING  SEQUENCE  AND  DATA  ACQUISITION 

The  Pulsar  50Q  pulser  was  placed  within  a  few  (<  10)  feet  from  the 
test  chamber  to  eliminate  the  difficulties  in  matching  the  load  of  the  cham¬ 
ber  to  a  long  coaxial  cable. 

The  throughput  delay  from  command  to  firing  of  the  pulser  is  250- 
350  ns,  with  ±  2  ns  jitter.  The  delay  of  AURORA  from  the  Marx  erection 
trigger  is  1.3  us.  We  were  told  this  delay  had  a  ±  50  nsec  jitter  but  (for 
constant  charging  parameters)  we  observed  far  less.  The  width  of  the  AURORA 
pulse  is  -  120  ns  at  FWHM.  To  insure  that  the  pulser  is  on  during  AURORA, 
the  pulser  was  fired  "just"  prior  to  the  AURORA  pulse  (i.e.,  «5  ns). 

The  pulser  was  fired  using  the  timing  circuit  in  Figure  13.  The 
position  of  the  pulser  firing  relative  to  the  AURORA  radiation  pulse  can  be 
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adjusted  using  delay  generator  I.  Delay  generator  II  is  used  to  simultane¬ 
ously  trigger  all  the  oscilloscopes  prior  to  both  the  pulser  and  AURORA 
f i ring. 

The  data  was  acquired  with  «25  oscilloscope  channels  and  photo¬ 
graphs  of  the  NIL.  Some  signals  were  redundantly  recorded  at  different 
gains  and/or  sweep  speeds. 

The  pulser  and  instrumentation  grounding  was  performed  using  the 
circuit  in  Figure  14.  The  cable  carrying  the  trigger  signal  to  the  10  A  was 
placed  outside  the  overall  shield  and  connected  to  the  trigger  output  in  the 
data  room  through  a  pulse  transformer.  This  was  done  because  a  high  voltage 
spike  was  discovered  on  that  cable  when  the  50  Q,  100  kV  pulser  fired  (i.e., 
the  Pulsar  50Q  sent  a  noise  pulse  back  down  its  own  trigger  line  to  the 
instrumentation  room  =<300  ns  after  it  received  a  trigger  from  the  Pulsar 
10A) .  When  the  cable  was  in  the  shield  it  injected  noise  on  the  signal 
cables. 
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SECTION  3 

EXPERIMENTAL  RESULTS 


The  experimental  results  will  be  presented  in  two  sections:  diag¬ 
nostics  and  nuclear  lightning  results. 

3.1  DIAGNOSTICS 

3.1.1  100  kV  Pulser  Behavior 

Scope  photos  of  the  voltage  on  the  pulser  with  and  without  the 
AURORA  radiation  pulse  are  shown  in  Figures  15A  and  15B.  The  pulser  oscil¬ 
lates  because  it  is  not  terminated  in  its  characteristic  impedance  (50  ft). 

It  is  initially  terminated  \in  essentially  an  open  circuit.  In  Figure  15B 
the  voltage  is  clamped  at  around  35  kV  during  the  AURORA  radiation  pulse. 

It  appears  that  a  large  current  is  flowing  somewhere  in  the  circuit.  The 
pulser  has  an  internal  inductance  of  1  uH.  The  internal  reactance  and  the 
external  resistance  to  ground  act  as  a  voltage  divider.  Exactly  where  the 
current  is  flowing  to  ground  is  unknown.  The  currents  that  are  observed  on 
the  current  ground  plates  are  not  sufficient  to  cause  such  a  voltage  drop. 
Conceivably  the  current  flowed  somewhere  within  the  pulser.  It  is  clearly 
not  an  arc  within  the  pulser,  because  an  arc  would  continue  to  discharge  the 
pulser  after  the  AURORA  pulse  ceased.  Our  best  explanation  is  to  assume 
that  the  current  flowed  through  dosed  conductive  air.  The  troube  with  this 
explanation  is  that  the  pulser,  to  the  best  of  our  knowledge,  was  completely 
filled  with  SFe  which  shouldn't  be  conductive  under  these  conditions.  Only 
if  it  were  not  completely  filled,  could  this  be  an  explanation  for  the 
observation.  As  indicated  in  Section  2  we  expected  the  pulser  to  generate 
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100  kV.  To  compensate  for  the  lower  voltage,  we  adjusted  the  delay  of  the 
pulser  so  that  its  natural  oscillations  would  place  it  on  a  voltage  peak 
during  the  radiation  pulse.  This  generated  the  largest  possible  voltges. 

The  final  timing  was  such  that  the  pulser  fired  =5  ns  prior  to  the  radiation 
pulse.  Notice  that  such  precise  timing  is  only  possible  if  AURORA'S  operat¬ 
ing  parameters  are  the  same  for  every  shot. 

3.1.2  Current  Probe  Response 

During  the  tests  of  the  current  plates,  we  discovered  that  they 
were  inductively  coupled  to  one  another.  If  current  is  driven  onto  the  in¬ 
ner  plate,  voltages  will  show  up  on  both  the  inner  plate  (pin)  and  the  outer 
plate  (annular).  Those  voltages  are  approximately  equal.  If  a  current  is 
driven  onto  the  outer  plate,  the  voltage  on  the  inner  plate  will  be  1/100  of 
that  observed  on  the  outer  plate.  This  means  that  to  determine  the  current 
on  the  outer  plate  the  current  on  the  pin  or  inner  plate  must  be  subtracted 
from  that  observed  on  th  outer  plate.  The  current  measured  on  the  inner 
plate  (pin)  is  virtually  unaffected  by  currents  on  the  outer  plate.  The 
response  of  the  current  plates  during  a  shot  in  which  nuclear  lightning  was 
observed  is  shown  in  Figure  16A  and  16B.  An  oscillation  with  a  frequency  of 
=•40  MHz  is  observed  on  these  signals.  We  attribute  this  to  LRC  oscillations 
of  the  current  plates. 

3.1.3  Grid  Response 

In  Figure  17  the  response  of  the  grid  and  a  voltage  probe  on  the 
pulser  output  are  shown.  The  voltage  indicated  on  the  grid  is  much  smaller 
than  expected.  The  voltage  observed  on  the  grid  is  apparently  due  to  capa¬ 
citive  coupling  to  the  cathode.  It  is  not  known,  at  this  time,  why  the  grid 
was  not  charged  up  by  the  large  current  («1  kA)  passing  through  it. 


Upper  trace: 

pin  current  111  A/d 
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Figure  16A.  Pin  current  with  AURORA  and  100  kV  pulser. 
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Figure  16B.  Large  ground  plate  current. 
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Figure  17.  Grid  Response 
3.1.4  AURORA  Diagnostics 

The  behavior  of  the  AURORA  radiation  pulse  was  monitored  using 
TLD's  and  gamma-dot  sensors  provided  by  the  AURORA  facility.  TLD's  were 
placed  in  10  to  15  locations  on  the  test  fixture  as  shown  in  Figure  18.  The 
total  dose  and  peak  dose  rate  determined  from  the  TLD  data  for  several  shots 
are  shown  in  Table  2.  The  dose  rate  of  a  typical  AURORA  shot  is  shown  in 
Figure  19.  The  TLD's  placed  inside  the  cathode  are  closest  to  the  location 
of  the  arc.  They  indicate  a  peak  dose  rate  of  9-10  x  10l°  rads/s. 


Figure  19.  AURORA  Pulse  Shape 


3.1.5  Noise  Background  Shots 

Several  background  shots  were  performed  to  determine  noise  levels 
on  the  signals.  These  include  shots  where  AURORA  was  fired  without  the  pul  - 
ser  and  shots  where  the  pulser  fired  without  AURORA.  In  Figure  20  the 
response  of  a  photodiode  to  the  AURORA  radiation  pulse  is  shown.  This 
response  has  been  attributed  to  fluorescence  of  the  plexiglass  structure  or 
air.  Photographs  taken  of  the  plexiglass  structure  during  the  AURORA  radia¬ 
tion  pulse  indicate  a  bright  fluorescence.  This  APO  signal  should  be  com¬ 
pared  with  HDL's  pulse  shape  monitor  of  Figure  19  or  the  cable  photocurrent 
shown  in  Figure  15B.  In  Figure  21  a  photograph  of  the  region  where  the  NIL 
will  occur  is  shown.  The  blurry  glow  is  attributed  to  the  plexiglass  struc¬ 
ture.  The  optical  system  was  designed  to  relay  the  image  of  the  1  cm  long 
NIL.  However,  an  image  with  a  smaller  magnification  will  be  relayed  by  the 
system  if  the  parabolic  mirror  is  ignored.  Although  thisimage  is  out  of 
focus,  it  is  sufficient  to  generate  the  blurry  background  on  the  photograph 
and  the  observed  photodiode  response.  The  signals  observed  on  the  current 
and  voltage  probes  during  shots  where  nuclear  lightning  was  observed  are 
several  volts  in  amplitude.  The  noise  signals  observed  on  both  AURORA  and 
100  kV  background  test  shots  are  insignificant  when  compared  to  the  observed 
signal  levels. 


Upper  trace:  extr 

Lower  trace: 

APD#1  50  mV/div 


100  ns/div 

Figure  20.  APD  response  to  AURORA  without  HV  pulser. 


Figure  21.  Plexiglass  structure  fluorescence  as  viewed 
by  the  optical  system. 
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3.1.6  Pre-Test  Shots  Without  AURORA 

Table  3  is  an  annotated  index  of  all  test  activities  including 
several  shots  without  the  AURORA  radiation  pulse  but  with  the  100  kV  pul- 
ser.  The  cathode  was  placed  close  enough  to  the  pin  for  arc-over  to  occur. 
When  arc-over  occured  the  pulser  dumped  all  of  its  energy  through  the  arc. 
This  is  shown  in  the  current  and  voltage  records  in  Figure  22.  This  arc  was 
much  brighter  than  the  NIL  event  since  it  draws  a  larger  current.  Figure  23, 
the  photodiode  response,  shows  the  photodiodes  saturating  immediately  after 
the  arc  closes  the  gap.  A  photograph  of  the  arc  is  shown  in  Figure  24. 

This  photograph  was  taken  with  the  aperture  set  at  its  smallest  value.  The 
wide  band  of  light  that  appears  to  be  in  the  background  is  an  image  of  the 
arc.  An  image  of  the  cathode  grid  and  pin  is  also  observed.  This  is  due  to 
the  accidental  (pin  hole)  imaging  properties  of  the  optical  system  discussed 
in  Section  3.1.5.  The  cathode  and  pin  apppear  as  a  shadow  in  front  of  the 
arc  because  of  the  small  aperture.  These  test  shots  verified  that  the 
optical  system  and  the  related  electronics  were  operating  properly. 

3.2  NUCLEAR  LIGHTNING  RESULTS 

To  accomodate  for  the  lower  voltage  available  from  the  pulser  dur¬ 
ing  the  AURORA  radiation  pulse,  the  separation  between  the  cathode  and  the 
pin  was  reduced.  This  increased  the  electric  field  to  encourage  the  crea¬ 
tion  of  a  nuclear  discharge.  Evidence  of  nuclear  lightning  was  first  ob¬ 
served  at  a  pin-cathode  spacing  of  10  cm  (5  shots).  Arc-over  occurred  at  a 
5  cm  spacing.  An  arc  «0.5  cm  long  was  observed  at  a  spacing  of  7.5  cm  (two 
shots).  Open  shutter  photographs  of  nuclear  streamers  observed  at  this 
spacing  are  shown  in  Figure  25A  and  B.  While  it  is  difficult  to  determine 
the  diameter  of  this  streamer  it  is  clearly  no  greater  than  1  mm.  The 
response  of  four  APD's  to  a  nuclear  streamer  in  shot  3715  (Figure  25B)  is 
shown  in  Figure  26.  APD  number  4  is  focused  on  the  region  at  the  tip  of  the 
pin.  APD  number  3  is  focused  on  the  next  region,  with  APDs  number  2  and  1 
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Annotated  table  of  experimental  activities  for  the  AURORA  nuclear  lightning  test 

(set-up  experiment  3/20-22/82) 
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Upper  trace:  (enhanced) 
voltage  78  kV/div 

Lower  trace:  extra 


200  ns/div 


Upper  trace:  extra 

Lower  trace: 

current  5.5  kA/div 


100  ns/div 


Figure  22.  Voltage  and  current  behavior  with  arc-over. 


Upper  trace:  extra 

Lower  trace: 

APD#1  500  mV/di v 


100  ns 


Upper  trace: 

APD#2  500  mV/di v 

Lower  trace: 

APD#3  500  mV/di v 


100  ns 

Figure  23.  APD  response  to  arc-over. 


“O  -o 

QJ  CD 
u  u 
c  c 

T3  03 
-£T  JZ 
C  C 

o>  o> 


>  > 

•r*  •  r— 

T3  -O 


3> 

-O 


o 

uO 


fT3  r— • 

S_ 

+->  O 
X  CL  > 
QJ  <t 

■O 


CL)  CD 

o  a 

rO  03 

s-  s- 


to 

c 


o 

o 


a»  qj 
a-  3 

Q-  O 


« — 

u 

“D 

L- 

CD 

03 

U 

c 

S- 

03 

03 

-C 

CD 

c 

r— 

CD 

a 

c 

> 

•r— 

o 

"O 

4-> 

> 

CD 

E 

to 

c 

o 

o 

o 

Cl. 

t-H 

to 

> 

QJ 

*d- 

•r— 

S-. 

=n= 

■O 

Q 

Q 

CL 

to 

CL 

C 

c 

C 

o 

•  • 

o 

• 

CD 

r-H 

to 

O 

CNJ 

03 

i- 

QJ 

4-> 

S- 

13 

S~ 

cn 

CD 

•r— 

Q. 

Li_ 

Q- 

46 


focused  on  the  following  two  regions.  The  field-of-view  of  each  region  is 
0.25  cm  long.  The  nuclear  streamer  appeared  about  150  ns  after  onset  of 
radiation.  When  the  image  of  the  NIL  enters  a  region  it  is  indicated  by  the 
signal  departing  from  the  bell  shaped  radiation  signal.  The  streamer  was 
observed  by  diodes  2-4,  but  never  entered  the  field  of  view  of  diode  1. 

The  onset  of  the  radiation  pulse  (indicated  by  the  fluorescence)  was  used  as 
a  fiducial  mark  to  correlate  the  time  on  the  traces. 

The  position  of  the  tip  of  the  NIL  is  plotted  as  a  functon  of  time 
in  Figure  27.  Timing  uncertainty  is  primarily  related  to  our  ability  to 
distinguish  the  onset  of  the  NIL  signal  in  the  scope  photograph  and  is  esti¬ 
mated  to  be  »±  6  ns.  The  location  uncertainty  of  the  photodetector  field  of 
view  with  respect  to  the  tip  of  the  pin  is  estimated  to  be  ±  0.13  cm.  The 
field-of-view  separations  have  no  uncertainty  since  all  detectors  share  a 
common  optical  path.  In  this  investigation  no  streamer  was  ever  long  enough 
to  enter  the  field-of-view  of  avalanche  photodiode  1.  The  peak  optical 
intensity  typically  occurred  175-260  ns  after  onset  of  radiation  and  voltage 
(or  «25-110  ns  after  initiation  of  the  streamer).  The  peak  intensity  ap¬ 
proximately  coincided  with  the  cessation  of  radiation.  That  the  streamer 
existed  at  these  late  times  but  was  not  seen  by  APD  1  is  clear  evidence  that 
it  stopped  propagating  before  entering  the  field  of  view  of  detector  1. 

This  observation  helps  confirm  these  were  nuclear  streamers  which  require 
conductive  air  rather  than  electrostatic  air  breakdown  "arcs". 

In  Figure  28  the  dose  rate,  the  applied  voltage,  the  pin  current, - 
and  the  response  of  APD  number  4  are  plotted  on  the  same  time  base  (high 
frequency  (»1  GHz)  noise  has  been  filtered  out  of  the  APD  data).  Four  times 
are  highlighted.  The  application  of  high  voltage  was  synchronized  with  the 
onset  of  radiation  (t0).  Notice  the  applied  voltage  dropped  as  dose  rate 
increased  and  recovered  on  cessation  of  radiation  .  The  peak  in  the  "fluo¬ 
rescence"  signal  coincides  with  the  peak  dose  rate  (t2).  The  nuclear  light¬ 
ning  streamer  enters  the  APD  field-of-view,  or  its  intensity  equals  the 
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fluorescent  signal,  at  tq  and  reaches  an  intensity  maximum  at  t  .  Notice 
there  is  apparently  an  incubation  time  for  the  formation  of  a  streamer. 

Dose  rate  and  current  are  actually  falling  as  the  streamer  intensity  grows. 
Applied  voltage  grows  as  dose  rate  and  current  decrease.  The  peak  intensity 
is  thought  to  occur  when  the  heating  and  cooling  rates  are  equal.  The  cur¬ 
rent  flowing  through  the  streamer  at  this  time  is  the  self-sustaining  cur¬ 
rent.  The  self-sustaining  current  appears  to  be  less  than  100  Amps.  The 
slow  decay  of  the  intensity  is  related  to  the  cooling  rate.  The  applied 
voltage  is  increasing  du'ing  the  growth  of  the  NIL.  The  diagnostics  of  this 
shot  are  very  similar  to  those  for  other  shots.  The  peak  current  on  shots 
where  nuclear  streamer  were  observed  ranges  from  100  A  to  440  A.  The 
applied  voltage  at  the  peak  of  AURORA  varied  from  «  24  kV  to  «  48  kV  between 
the  shots. 

3.3  COMPARISON  OF  EXPERIMENTAL  RESULTS  WITH  THEORY 

The  fact  that  dose  rate  current  and  voltage  were  all  varying  si¬ 
multaneously  during  the  streamer  propagation  make  it  difficult  to  compare 
the  instantaneous  velocity  with  expectations  based  on  theory.  The  line 
drawn  through  the  data  points  of  Figure  27  is  just  for  comparison  and  does 
not  represent  a  true  "fit"  to  the  data. 

In  Reference  2,  Frese,  et  al.,  did  a  preliminary  estimate  of  the 
propagation  theory  of  nucler  lightning.  In  this  oversimplified  model  they 
assumed  the  electric  field  was  a  constant  Ea  (the  avalanche  field)  within 
a  spherical  shell  of  radius  (R),  so  that  the  current  density  increased  with 
decreasing  radius  (r)  as  r-2.  They  further  assumed  the  power  P  =  /JEadA, 
deposited  in  a  spherical  volume  of  t  '.<lius  (a)  would  adiabatical  ly*  heat  the 
air  contained  therein.  Calculations  suggest  only  radiation  losses  can  occur 
on  these  time  scales.  Expansion  and  conduction  don’t  have  time  to  move  much 
energy.  Assuming  a  constant  channel  diameter,  which  seems  consistent  both 

*  without  energy  loss  mechanisms 
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with  theory  and  observation  for  this  short  length,  the  streamer  will  grow 
longitudinally  at  a  rate  determined  by  the  ratio  of  power  deposited  to  heat 
(Q)  required  to  thermally  ionize  the  air.  According  to  this  model  the 
velocity  v  is 


v  =  o0Ea2R2/Qa  (1) 

where  aQ  is  the  radiation  induced  conductivity  of  air,  aQ=>KD.  Frese,  et  al., 
assumed  Q  =  40  J/cm2  and  deduced  a  =  0.025  cm  by  taking  v  =  1  x  107  cm/s 
from  the  Mike  Shot  NIL  observations. 

In  this  experiment  we  have  direct  measurements  of  the  applied 
voltage,  current  and  streamer  displacement.  Dose  rate  is  easily  and  accu¬ 
rately  derived  from  dosimetry.  Air  conductivity  can  be  approximated  from 
dose  rate  and  published  data.  Unfortunately  all  these  parameters  vary  rap¬ 
idly  over  the  50  ns  that  the  observed  streamer  propagated.  However  the  dis¬ 
placement  X  of  the  channel  tip  can  be  estimated  and  compared  with  expecta¬ 
tions.  Alternatively  we  can  treat  Qa  as  an  unknown  and  use  these  experimen¬ 
tal  data  to  derive  Qa. 

First,  the  total  instantaneous  power  dissipated  in  the  air  cell  is 
V(t)I(t).  Solution  of  the  electrostatic  problem  where  E  =  £  out  to  radius 

a 

R  and  falls  as  r-2  there  after  results  in  V(R)  =  Vapp-|/2  (Ref.  4). 

Therefore  a  channel  of  constant  diameter  can  propagate  no  farther 

than 


1 

2ira2Q 


t 

/  V(t)  I (t )  dt 
o 


X(t)  = 


(2) 


This  is  certainly  an  overestimate  because: 

1)  The  adiabatic  approximation  ignores  all  power  loss  mechanisms 
(i.e.,  radiation). 

2)  Most  (90%)  of  the  power  actually  flows  through  a  cylinder 
surrounding  the  pin  and/or  channel,  not  just  across  a  hemi¬ 
sphere  surrounding  the  tip. 

Another  way  to  estimate  the  channel  propagation  from  the  available  data 
(V(t)  and  D(t) )  is  by  substituting  for  Ea  and  R  in  terms  of  the  electrode 
potential  of  the  experiment  V,  (the  voltage  at  the  surface  of  the  avalanche 
region  is  V/2)  (Ref.  4).  For  this  case  the  position  as  a  function  of  time 
is  given  by 

t  . 

X(t)  =  JL  /  0  V2  dt  (3) 

4Qa  o 

where  K  is  the  constant  relating  dose  rate  D  and  air  conductivity  a, 

(K  «Kq  e  u/a  =  7  x  10-15  (n  cm  rad/s))-1. 

Figure  29  shows  X(t)  calculated  from  Equation  2  and  3  as  well  as 
the  measured  data  for  event  #3715.  For  comparison  purposes  the  data  is 
plotted  as  it  would  appear  if  Q  =  40  J/cm3  and  a  =  0.025  cm  (values  taken 
from  Reference  2).  It  would  appear  that  these  data  can  be  accomodated  by 
adiabatic  theory  if  Qa  is  considerably  less  than  inferred  from  previous  nuc¬ 
lear  tests  (i.e,  Qa  about  .05  J/cm2  rather  than  about  1  J/cm2).  We  suppose 
that  the  channel  diameter  2 a  of  a  0.75  cm  long  streamer  may  be  considerable 
smaller  than  the  diameter  of  a  6  x  IQ4  cm  long  streamer.  The  plot  of  Equa¬ 
tion  2  indicates  that  the  measured  data  is  physically  reasonable  from  an 
energy  viewpoint. 
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The  reason  for  the  apparent  ( » 125-1 50  ns)  incubation  time  sug¬ 
gested  by  these  data  is  not  known.  It  is  conceivable  that  at  early  times 
when  the  tip  is  near  the  metal  pin  that  the  energy  losses  are  larger  than 
after  the  tip  has  moved  away.  Thus  the  Qa  product  is  so  initially  large  and 
hence  the  tip  velocity  so  low  that  the  tip  remains  under  the  cooling  influ¬ 
ence  of  the  metal  pin  for  roughly  125  ns  (the  tip  velocity  should  be  propor¬ 
tional  to  the  di fference  between  the  heating  and  cooling  rates). 
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SECTION  4 
SUMMARY 


This  report  describes  an  experiment  designed  to  duplicate  the  con¬ 
ditions  at  the  tip  of  nuclear  lightning  in  a  laboratory  environment  by  scal¬ 
ing  dose  rate  up  and  dimensions  down.  All  indications  are  that  the  experi¬ 
ment  produced  miniature  nuclear  lightning  bolts  (about  0.75  cm  long)  at 
least  5  times  in  an  AURORA  exposure.  These  events  were  photographed,  elec¬ 
trically  characterized  and  the  discharge  propagation  was  measured.  Agree¬ 
ment  with  expectations  based  on  previous  atmospheric  test  observations  and 
theory  is  judged  to  be  reasonable.  This  experiment  lends  considerable  cre¬ 
dibility  to  the  propagation  theory  of  nuclear  induced  lightning. 

There  remain  some  uncertainties,  in  the  theory  of  nuclear  lightn¬ 
ing  particularly  in  the  general  areas  of  air  chemistry  (i.e.,  species  and 
reaction  rates)  and  radiative  transfer.  That  nuclear  lightning  can  be 
created  in  laboratory  conditions  has  been  demonstrated.  Experiments  such  as 
time  resolved  spectroscopy  designed  to  reduce  the  remaining  questions  should 
be  considered. 
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APPENDIX 

ELECTRIC  FIELD  MEASUREMENT  (LANGMUIRE  PROBE) 


The  boundary  layer  voltage  is  measured  by  placing  a  grid  «  2  cm  below 
the  cathode.  The  grid  is  made  of  4  mil  wire  spaced  1  cm  apart.  The  grid 
can  best  be  described  in  this  case  as  a  Langmuire  probe.  As  electrons  pass 
through  the  grid,  some  electrons  will  strike  the  grid  and  charge  it  up  until 
the  voltage  is  slightly  larger  than  the  ambient  electric  potential,  then  the 
incident  electrons  will  be  deflected  from  the  grid.  Loeb  (Ref.  Al)  shows 
that  this  votage  increase  is  given  by 

v  =J^m  I=5± 

W  2e  T+M. 


Where  M+  is  the  atomic  mass  of  positive  ions,  M_  is  the  atomic  mass  of  the 
electrons,  T_  and  T+  are  the  temperature  of  the  electrons  and  the  ions,  re¬ 
spectively.  k  is  boltzmann's  constant  and  e  is  the  electro  charge.  T_  for 
“  1  eV  electrons  is  11.6  x  103°  K,  T+  for  N2-02  ions  is  -  300°  K,  M_  is 
5.5  x  10-4  AMU  for  electrons,  M+  for  N2  and  /or  02  is  «  60  AMU,  e  is  1.6  x 
10-19  coul  and  k  is  1.38  x  10~23  J/mole°  K.  For  this  case  Vy  is  “8  volts. 
This  is  insignificant  relative  to  the  90  kV  being  measured. 

The  entire  grid  is  10  cm  in  diameter.  The  percent  of  the  current  avail¬ 
able  to  the  grid  is  proportional  to  the  area  subtended  by  the  grid,  shown  by 


%  of  current  on  grid  = 


cross  sectional  area  of  grid  wire 
total  area  subtended  by  the  grid 


x 


100 


Al.  L.B.  Loeb,  Basic  Processes  of  Gaseous  Electronics,  University  of  Cali¬ 
fornia  Press,  Berkeley,  California,  1955,  p.337. 
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For  the  grid  described  above,  the  area  of  the  wire  is  1.72  x  10-4  m2,  and 
the  total  area  is  7.85  x  10"3  m2,  resulting  in  2%  of  the  total  current 
picked  up  by  the  grid.  Most  of  the  current  flowing  to  the  pin  will  also 
flow  through  the  grid. 

Currents  ranging  fom  100  to  1,000  Amps  were  expected  at  the  pin.  Grid 
currents  between  2  and  20  amps  were  expected.  The  voltage  probe  used  to 
monitor  the  grid  could  not  draw  more  than  the  grid  current,  or  the  grid 
voltage  would  drop  below  the  value  of  the  ambient  field  potential  near  the 
grid.  The  voltageon  the  grid  was  measured  with  a  -  28  kft  voltage  divider. 
The  signal  was  carried  to  the  divider  via  an  insulated  cable  placed  approxi¬ 
mately  along  an  equipotential .  The  voltage  divider  consisted  of  20,  1.4  kft 
resistors  and  one  1.8  ft  measurement  resistor.  The  response  of  the  divider 
was  calculated  from 


w  -  w  L 

V0UT  ~  "IN  rl  +  rd 

where  Rl  (load  resistor)  is  1.8  ft  and  Rq  (divider  resistor)  is  28  kft. 
Vqut/Vin  =  x  10-5.  For  a  grid  voltage  of  a  90  kV,  the  output  volt¬ 
age  5.8  V  and  the  current  drawn  amps  response  time  of  this  probe  was  given 
by  (Ref.  5) 


TP 


1,2  RS  C5 


where  N  (number  of  resistors)  was  20,  Rs  (resistor  value)  was  1.4  kft  and 
Cs  (shunt  capacitance  per  resistor)  was  *  0.5  pF.  For  this  probe  tp  = 

28  ns.  The  voltage  across  the  1.8  ft  resistor  was  sent  via  an  RG-58  (50  ft) 
coaxial  cable  to  the  data  room  .  The  grid  and  voltage  probe  are  shown  in 
Figure  Al. 
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RG-58  cable 


sensor. 
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